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摘 要 :确定 性 横向 迁移 (deterministic lateral displacement,DLD ) 是 一 种 常用 的 微粒 分 离 技术 ,依据 

粒子 尺寸 和 微 流 道 层 流 特征 分 高 颗粒 。 传 统 的 DLD 理论 一 般 只 考虑 刚性 粒子 ,难以 描述 细胞 等 柔 

od 离 过 程 。 结 合 分 子 动力 学 粗 粒度 模型 和 雷诺 相似 性 原理 ,对 DLD 装置 中 柔性 细胞 的 分 
形 和 受 力 情况 进行 了 研究 。 仿真 结 果 发 现 , 纳 观 尺度 下 的 DLD 细胞 分 离 路 线 与 宏观 尺度 的 

离 原理 相符 。 而 在 微观 尺度 下 ,细胞 分 离 过 程 中 细胞 的 变形 破坏 主要 出 现在 细胞 与 微 柱 接 

触 处 ,而 且 该 处 所 受 应 力 最 大 。 研 究 结 果 揭 示 了 柔性 颗粒 通过 DLD 装置 的 运动 细节 。 
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Molecular dynamics study of cell motion characteristics in 
deterministic lateral displacement devices 
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Abstract ; Deterministic Lateral Displacement ( DLD) is a common particle separation technique in which 
the particles are separated according to their size and the characteristics of microfluidic laminar flow. The 
traditional DLD theory only considers rigid particles and it is difficult to describe the separation process of 
flexible particles such as cells. In this paper ,the separation , deformation and stress of flexible cell in DLD 
device were studied by using coarser particle molecular dynamics model and Reynolds similarity principle. 
The simulation results show that the separation route of DLD cells at nano scale is consistent with the 
general separation principle at macro-scale. However, in the microscopic scale, the deformation and 
destruction of cells in the process of cell separation mainly occurred at the contact point between cells and 


microcolumn , where the maximum stress was exerted. The results reveal the motion details of the flexible 
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近年 来 , 微 流 控 技 术 发 展 迅速 ,已 成 为 许多 工程 
和 生物 医学 应 用 的 重要 组 成 部 分 站。 微细 颗粒 分 
选 是 该 技术 的 关键 应 用 之 一 。 目 前 常用 的 高 通 量 微 
粒 分 选 方法 主要 分 为 主动 分 选 法 和 被 动 分 选 法 。 被 
动 分 选 法 利用 粒子 , 流 场 和 通道 结构 之 间 的 相互 作 
用 实现 ,不 需要 外 部 场 *“。 确 定性 横向 迁移 ( deter- 
ministic lateral displacement, DLD ) 是 其 中 一 种 被 动 
分 选 方法 , 因 其 操作 简单 .时 间 短 、 分 离 的 微粒 范围 
广 及 能 够 实现 连续 分 离 等 优点 而 被 广泛 研究 5 。 
DLD 分 离 技 术 是 指 依据 粒子 尺寸 ,利用 其 微 通 道内 
特定 设置 的 微 柱 精 确 控 制 微粒 的 轨迹 ,促进 不 同 大 
小 粒子 分 离 的 技术 [9 。 传 统 的 DLD 理论 一 般 只 
虑 刚性 粒子 ,然而 实际 应 用 中 ,DLD 常常 用 来 分 离 生 
物 颗粒 ,比如 细胞 等 ,具有 很 强 的 柔性 ,也 存在 多 种 
形状 ,传统 DLD 分 选 理论 很 难说 明 这 类 问题 。 

目前 国内 外 对 细胞 DLD 的 分 离 过 程 已 经 开展 
TERT, BEECH 55" 依据 剪 切 率 和 变形 
性 改变 细胞 的 有 效 尺 寸 ,从 而 确定 分 离 相关 参数 。 
HOLM 等 :根据 分 离 过 程 中 细胞 不 同 取向 导致 的 
不 同 有 效 尺寸 来 改变 分 离 装置 ,从 而 实现 细胞 分 离 。 
以 上 研究 局 限于 分 析 细 胞 尺寸 和 变形 性 对 分 离 过 程 
的 统计 性 影响 ,对 于 具体 细胞 在 流 场 中 的 运动 变形 
及 其 力学 特性 的 实验 研究 , 受 限 于 观测 手段 的 困难 
性 ,目前 相关 的 报道 仍然 比较 稀少 。 

数值 模拟 方法 在 柔性 细胞 的 DLD 分 选 研究 中 
发 挥 了 重要 作用 。KHODAEE 等 ("1 运用 有 限 元 法 研 
究 了 在 DLD 微 流 控 装置 中 循环 肿瘤 细胞 和 白细胞 
的 运动 . 受 力 和 分 离 , WEIL 等 "运用 浸入 边界 -格子 
玻 尔 效 曼 法 对 DLD 分 离 柔 性 细胞 进行 了 模拟 ， 
ZHANG 45 ^5 通过 耗 散 粒子 动力 学 方法 探究 了 微 
柱 形状 与 细胞 变形 之 间 的 作用 关系 。 以 上 研究 从 宏 
观 角 度 分 析 细 胞 的 变形 ,将 细胞 作为 一 个 柔性 整体 
进行 处 理 。 实 际 上 细胞 包括 细胞 膜 和 细胞 质 等 微观 
结构 ,这 些 结构 对 细胞 变形 和 力学 特征 的 影响 不 可 
忽视 ,这 方面 的 研究 非常 缺乏 。 分 子 动力 学 是 进行 
微观 模拟 的 有 效 方法 ,目前 已 经 存在 一 些 关 于 细胞 
膜 结构 和 组 成 的 分 子 动力 学 研究 ”|。BARNS 
等 2 结合 原子 力 显 微 镜 压 痕 和 粗 粒 度 颗粒 法 研究 
了 红细胞 的 力学 特性 , CHANG 等 (中 运用 MD/DPD 
多 尺度 框架 预测 了 健康 和 疾病 中 的 红细胞 形态 和 应 


力 。 受 限于 计算 效率 ,这 些 研究 一 般 只 关注 细胞 局 
部 的 结构 和 组 织 , 对 于 DLD 装置 中 整个 细胞 分 离 过 
程 的 分 子 动力 学 仿真 ,尚未 发 现 相关 报道 。 

鉴于 此 ,本 研究 提出 结合 粗 粒 度 分 子 动力 学 模 
型 和 雷诺 相似 性 原理 ,有 望 解决 DLD 中 柔性 细胞 运 
动 仿真 的 难题 。 根 据 忌 片 结 构建 立 二 维 仿真 模型 ， 
选取 粗 粒 化 磷脂 分 子 模型 ,编写 程序 对 具有 双 分 子 
层 的 细胞 膜 系统 进行 建 模 。 通 过 仿真 结果 分 析 了 各 
个 阶段 和 不 同 微 柱 直径 情况 下 ,细胞 在 DLD 芯片 中 
的 运动 . 受 力 和 变形 机 制 。 


1 仿真 模型 及 参数 


1.1 DLD 芯片 模型 及 参数 


典型 的 DLD 芯片 模型 如 图 1 所 示 , 在 一 个 微 流 
道中 ,周期 性 地 排列 了 若干 微 柱 ,每 列 微 柱 与 前 一 列 
在 纵向 稍微 错开 。 阵 列 的 横向 周期 为 Y, 微 柱 直 径 
为 4, 微 柱 间距 为 a, 阵列 偏 移 率 A = 1/N, 偏 移 角 度 
为 9, 相 邻 列 的 纵向 偏 移 量 为 a xA。 流 体 在 流 道内 
流动 时 ,形成 倾斜 的 主流 道 和 横向 的 副 流 道 ,主流 道 
宽度 等 于 微 柱 间距 , 即 为 a, 副 流 道 宽度 为 5。 因 微 
柱 阻 碍 效果 ,主流 道上 的 流体 按 纵向 方向 被 分 成 若 
干 条 有 具有 相同 流量 的 流 层 ,其 临界 偏转 半径 为 R.， 
当 粒 子 直径 小 于 临界 偏转 直径 2R. 时 ,流体 中 的 粒 
子 会 穿 过 副 流 道 到 达 下 一 主流 道 某 一 流 层 ,然后 继 
续 沿 主流 道 方向 流动 ,呈现 ”ZZ 走向。 反之 , A 
子 直径 大 于 2R. 时 ,粒子 在 运动 过 程 中 ,每 经 过 一 个 
微 柱 后 都 会 迁移 到 下 一 个 主流 道 的 同一 类 流 层 中 ， 
呈现 “LD? 走 向 。 根 据 粒 子 直 径 尺 十 的 不 同 ,粒子 在 
DLD 芯片 中 随 流体 流动 到 相应 收集 通道 ,最终 实现 
不 同 大 小 粒子 的 分 离 。 

由 于 计算 量 的 限制 ,本 研究 根据 雷诺 相似 准则 ， 
采用 了 等 比例 缩小 的 仿真 模型 。 保 证 模型 内 的 流体 
与 实际 DLD 内 流体 的 雷诺 数 相同 , 即 

Patiala Punka -及 (1) 

从 Mm 
其 中 :pu us La 和 Hu 分 别 为 实际 DLD 芯片 中 流体 的 
密度 .速度 .特征 长 度 和 黏 性 系数 ;pw Um Lu 和 An 
分 别 为 分 子 动力 学 模型 中 流体 的 密度 .速度 .特征 长 
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度 和 条 性 系数 。 调 整 模 型 的 流 道 尺 寸 与 流速 等 各 项 
参数 ,使 得 模型 中 流体 与 实际 DLD 芯片 中 流体 的 雷 
诺 数 一 致 。 根 据 本 研究 设置 ,模型 的 雷诺 数 为 0.2 ~ 
0.3. 参照 文献 [6 ] ,为 获得 较 好 的 抛物 线 流 线 、 较 
为 分 明 的 细胞 分 离 特征 和 形变 特征 ,确定 微 柱 轮 廊 
为 圆柱 形 ,并 设 定 模型 各 项 参数 如 表 1 所 示 。 


图 1 DLD 芯片 模型 
Fig.1 DLD chip model 
R1 DLD 模型 参数 表 
Tab.1  DLD model parameters table 


DLD 芯片 /nm? N a/nm b/nm d/nm 0/(°) 


126 x42 3 24 26.3 14 «22 18.43 


1.2 ”细胞 模型 及 参数 


构建 好 DLD 分 子 动力 学 模型 之 后 ,在 进行 DLD 
分 离 细 胞 的 分 子 动力 学 模拟 之 前 ,还 需 在 模型 中 加 
人 具备 一 定 结构 和 尺寸 的 细胞 模型 。 细 胞 的 结构 可 
认为 主要 是 由 细胞 腊 和 被 细胞 膜 包 庄 住 的 细胞 基质 
组 成 。 本 研究 主要 涉及 细胞 膜 的 流动 特性 及 其 在 高 
剪 切 率 下 的 受 力 ,参考 文献 [17] ,考虑 在 细胞 膜 建 
模 方 面 主要 对 磷脂 双 分 子 层 进行 建 模 。 在 DLD 分 
选 过 程 中 ,细胞 承受 较 强 的 流体 动力 学 应 力 和 边界 
的 机 械 反 作用 力 , 相 比 之 下 ,在 较 短 的 时 间 内 可 以 忽 
略 物质 扩散 和 流体 秋 度 影响 , 且 目 前 关于 细胞 基质 
的 粗 粒 化 模型 非常 少 , 故 将 细胞 基质 简化 成 由 水 分 
子 构成 的 细胞 内 液 ,同时 为 避免 渗透 压 影 响 ,DLD 中 
细胞 外 液 和 流体 也 为 水 。 由 此 建立 由 细胞 膜 和 细胞 
内 液 组 成 的 二 维 细胞 模型 ,如 图 2 所 示 。 

本 研究 采用 粗 粒 化 DOPC 磷脂 模型 ”对 细胞 
膜 进 行 建 模 ,DOPC 模型 如 图 2(a) 所 示 , 该 模型 能 
体现 磷脂 分 子 的 亲 玻 水 性 , 较 好 地 表达 细胞 膜 的 流 
动 特性 。 模 型 包括 15 个 粗 粒 化 的 粒子 ,由 带 正 电 的 
胆 碱 . 带 负电 的 磷酸 盐 . 带 极 性 的 甘油 形成 的 磷脂 头 


ten) 1*8 


y y 
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部 和 两 条 碳 氢 键 尾 链 构 成 。 对 于 水 分 子 , 也 采用 粗 
粒 化 模型 ,如 图 2 (b) Brzs ,每 个 粒子 代表 一 个 水 分 
子 , 它 的 极 性 特性 由 座 在 内 部 的 偶 极 子 实现 。 


(a) (b) (c) 


Oo Q o. 
d Qu Ò 
A A EO 
9 9 QO o 
© 60 O O Gy 
© O60 Q Q Os 
© OO OW Òr 
6 60 O' o0 
(e) (f) (g) (h) 


图 2 ”细胞 模型 
Fig.2 The cell model 
分 子 动力 学 细胞 模型 仿真 运算 的 基础 是 具备 完 
整 的 分 子 结构 和 选择 合适 的 势 函 数 。 构 建 细胞 膜 结 
构 模 型 的 基本 做 法 是 : 先 根据 DOPC 模型 中 各 粒子 
的 数据 建立 单个 磷脂 分 子 ; 然 后 以 单个 磷脂 分 子 作 
为 复制 对 象 ,实现 磷脂 双 分 子 层 的 构建 ;之 后 依据 粗 
粒 化 水 分 子 模型 ,建立 图 2(c) 所 示 的 细胞 内 液 和 细 
胞 外 液 ,通过 以 上 步骤 实现 细胞 模型 的 建立 。 水 分 
子 和 磷脂 分 子 的 粗 粒 化 模型 如 图 2 所 示 , 图 中 圆 形 
粒子 中 的 数字 代表 原子 类 型 ,如 图 2(d)、(e) 所 示 ; 
粒子 与 粒子 间 的 短 横 线 代表 键 连接 , 键 旁 边 的 字母 
(a b ede) 代 表 键 类 型 ,如 图 2(f 所 示 ; 两 条 短 横 
线 之 间 的 角度 代表 键 角 ,希腊 字母 (a、B、y、5) 代 表 
键 角 类 型 ,如 图 2(g) 所 示 ; 通 过 复制 磷脂 分 子 以 及 
加 入 水 分 子 建立 分 子 动力 学 粗 粒 化 细胞 模型 ,如 图 
2(h) 所 示 。 
粗 粒 化 水 分 子 模型 和 DOPC 磷脂 模型 的 势 函数 
为 ORSI 等 .站 开发 的 ELBA 力 场 ,该 力 场 能 够 体现 
磷脂 模型 中 的 静电 、 偶 极 子 和 极 性 ,同时 能 够 体现 水 
分 子 的 极 性 。 粗 粒 化 水 分 子 间 的 作用 总 势能 ， 
U; =U; + UY (2) 


AA o 
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3 (u; : r;) (M; $ r;) 
a 5 ]s (4) 


其 中 ;U4 为 LI 8984 UP RV REA s e .eu vc 
分 别 为 势 井 深度 .真空 介 电 常 数 和 作用 势 为 0 时 两 


] Ai Ai 
Li = J 
2 vl r, 
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构成 一 个 粗 粒 化 细胞 ,黑色 圆 形 和 半圆 形 是 DLD is 
片 的 微 柱 , 细 胞 悬浮 在 左边 随 流体 从 左 至 右 流动 。 
本 研究 模拟 的 实际 物理 对 象 为 球形 红细胞 ,直径 约 
为 3 pm。 分 子 动力 学 模拟 需要 大 量 计算 资源 ,为 提 


粒子 间 的 距离 ;jw, yw 分 别 为 粒子 i 和 粒子 j 的 偶 极 
矩 ;r, 为 粒子 ; 中心 点 指向 粒子 ) 中心 点 的 向 量 ;s, 为 
由 转换 半径 和 截断 半径 确定 的 分 段 函数 。 

对 于 DOPC 磷脂 模型 ,其 粒子 间 总 势能 由 4 部 
分 组 成 :2 个 粒子 间 的 键 连接 能 .3 个 粒子 间 的 键 角 
能 .带电 粒子 间 的 静电 作用 能 和 带电 粒子 与 极 性 粒 
子 间 的 作用 能 ,分 别 为 


1 
U; s -lo)? (5) 
1 2 
Up = y cosa — cosa ) (6) 
Q,Q, ry 
00 . viv fa Py 
uU" = ma | (7) 
2 3 
Qu .. Q; e Fj Fy 
U; EL (2 a hn oO 


其 中 :kw NAER Z D; d oc on 分 别 为 实际 键 
长 参考 键 长 .实际 键 角 和 参考 键 角 ;0Q;、0; 分 别 为 粒 
Fi 和 粒子 7/ 的 电荷 量 。 

式 (2) ~ (8) 的 势能 函数 描述 了 粗 粒 度 模 型 中 
各 原子 间 的 相互 作用 形式 。 通 过 求解 势能 函数 的 梯 
度 , 可 以 得 到 数值 仿真 中 的 作用 力 和 加 速度 方程 。 
这 组 方程 确定 了 仿真 中 所 有 粗 粒 化 原子 的 微观 运动 
行为 ,宏观 统计 上 表现 为 本 研究 中 流体 的 流动 和 细 
胞 运动 。DOPC 磷脂 模型 势 画 数 中 的 各 项 参数 取 值 
参考 文献 [20 ] ,各 分 子 模型 的 各 项 参数 如 表 2 所 示 。 

R2 粗 粒 化 水 分 子 、 磷 脂 分 子 模型 参数 表 


Tab.2 Table of parameters of coarse-grained water molecule 


and phospholipid molecule model 


原子 名 称 偶 极 矩 / BE BIRRE 扩散 系数 / 


类 型 的 D C À (10-9m.s-2) 
1 水 分 子 0.541 0 3.0 2.7 

2 胆 碱 0 +0.7e 4.5 0.9 

3 磷酸 盐 0 -0.7e 5.2 1.2 

4 甘油 0.333 0 4.6 1.2 

5 酯 基 0.416 0 4.6 1.2 

6 碳 氧 键 0 0 4.5 0.9 


鉴于 DLD 结构 的 周期 性 ,本 研究 仅 选 一 个 周 
期 。 分 子 动力 学 模型 的 边界 处 , 均 施 加 周期 性 边界 


高 计算 效率 ,取得 有 效 计算 结果 ,本 研究 采用 了 相似 
性 原理 ,在 模拟 中 采用 了 缩小 约 100 倍 的 细胞 模型 。 
模型 细胞 直径 为 9 ~40 nm , 微 柱 直 径 为 14 ~22 nm, 
利用 大 规模 原子 /分 子 并 行 模拟 器 Lammps 进行 计 
算 ,然后 应 用 可 视 化 工具 ovito 对 模拟 结果 进行 可 视 
化 和 后 处 理 。 


JENN ——— 
图 3 DLD 芯片 的 分 子 动力 学 模型 
Fig.3 Molecular dynamics model of DLD chip 

模拟 的 具体 细节 如 下 :整体 采用 NVE 系 综 进 行 
模拟 ,使 用 恒温 模式 ,温度 控制 在 300 ,施加 外 势 场 
的 方法 驱动 流 场 流动 ,其 效果 对 应 于 宏观 流 场 下 的 
压力 梯度 。 用 光滑 的 固定 边界 模拟 微 柱 ,时 间 步 长 
为 0.1fs。 关 于 磷脂 模型 的 验证 在 参考 文献 [20 |] 中 
已 有 详细 描述 ,目前 Lammps 的 应 用 也 非常 成 熟 。 
本 研究 将 直接 根据 模拟 结果 进行 分 析 , 对 模型 和 算 
法 的 验证 不 再 著述 。 


2 结果 与 讨论 


2.1 细胞 运动 分 析 


按照 DLD 分 离 理 论 , 当 细胞 直径 大 于 临界 偏转 
直径 时 ,粒子 运动 模式 为 “LD” 模 式 , 而 细胞 直径 小 
于 临界 偏转 直径 时 , 则 呈现 “ZZ” 运动 模 式 。 在 本 模 
拟 中 , 当 微 柱 直 径 为 18 nm 时 ,由 理论 公式 IC 计算 得 
到 该 DLD 芯片 的 临界 偏转 直径 约 为 19. 8 nm ,细胞 
直径 为 27.6 nm 时 ,该 微 柱 直 径 下 模拟 结果 如 图 4 
(a) 所 示 , 细 胞 的 轨迹 为 “LD” 的 运动 模式 ; 当 细 胞 直 
径 为 9 nm 时 ,模拟 结果 如 图 4(e) 所 示 , 其 轨迹 呈现 
了 “ZZ” 走 向 ,两 种 模拟 结果 均 与 理论 结果 一 致 。 

将 细胞 在 微 流 道中 的 主要 行经 点 分 为 3 个 区 域 
( 工 工 亚 ), 图 4(b) ~(d) 分 别 为 对 应 细胞 运动 至 


条 件 PBC ,以 模拟 DLD 的 周期 性 结构 ,如 图 3 所 示 。 
模型 内 灰色 区 域 为 液态 水 , 左 侧 环形 及 其 内 部 的 水 


区 域 [、 工 ` 亚 时 的 情形 。 白 色 虚 线圈 标记 了 细胞 与 
微 柱 接触 时 的 情况 ,观察 发 现 ,在 接触 部 位 ,细胞 膜 
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与 微 柱 之 间 存 在 少量 流体 ,由 此 形成 一 层 液 膜 。 如 
图 4(p) 所 示 , 当 细胞 经 过 区 域 工 时 ,由 于 细胞 直径 
大 于 微 柱 间距 ,在 流体 的 推动 和 微 柱 挤 压 作 用 下 , 细 
胞 上 下 两 端 变形 ,但 由 于 液 膜 的 存在 , 微 柱 对 细胞 没 
有 明显 的 黏 滞 效果 。 当 细胞 到 达 区 域 开 时 ,如 网 4 
(c) 所 示 ,细胞 与 微 柱 碰撞 发 生 了 变形 ,靠近 微 柱 表 
面 的 细胞 膜 经 过 流体 的 挤 压 和 剪 切 作用 ,加 上 微 柱 
的 黏 滞 和 阻碍 作用 ,发生 了 显著 变形 。 同 时 ,在 微 柱 
的 阻碍 下 ,由 于 流体 惯性 ,以 及 细胞 膜 与 细胞 内 液 的 
刚度 差别 较 大 ,细胞 内 液 往 微 柱 一 侧 挤 压 , 细 胞 内 部 
出 现 空 从 。 这 是 模型 简化 引起 的 ,由 于 本 人 研究 将 细 
胞 质 简 化 为 水 ,忽略 了 秋 弹 性 ,忽略 了 细胞 质 与 细胞 
膜 之 间 的 结合 作用 ,因而 在 较 强 的 应 力作 用 下 ,有 可 
能 造成 细胞 膜 与 细胞 质 的 分 离 。 在 区 域 焉 处 ,如 图 
4(d) 所 示 , 细 胞 腊 也 与 微 柱 发 生 了 接触 ,由 于 其 速 
度 方向 已 经 发 生 了 较 大 变化 ,细胞 横向 迁移 ,此 时 微 
柱 主要 对 其 产生 么 清 作用 ,因此 细胞 表现 出 较 大 变 
形 , 细 胞 膜 有 轻微 破裂 。 


注 :a.b.c.d 细胞 直径 为 27.6 nm;e 细胞 直径 为 9 nm, 
图 4 细胞 在 DLD 芯片 中 的 运动 轨迹 与 结构 变化 
Fig.4 Cell trajectory and structure change in DLD chip 


根据 模拟 结果 中 各 分 子 的 空间 位 置 可 以 观察 细 
胞 的 结构 变化 ,通过 微 流 道内 流体 的 流 场 和 压力 分 
布 则 可 以 更 直观 地 理解 细胞 在 分 离 过程 中 的 运动 特 
性 。 如 图 5 所 示 , 图 5(a) ~5(c) 分 别 展示 了 细胞 经 
过 工 . 工 ` 亚 区 域 时 微 流 道 内 的 流 场 和 压力 情况 。 其 
中 ,细胞 的 实际 位 置 通过 黑色 线圈 表征 的 细胞 膜 展 
FR , 流 场 中 微 流 体 的 速度 通过 不 同 长 度 和 方向 的 白 
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色 矢 量 箭头 表示 , 微 流 道内 部 压力 则 是 根据 实时 压 
力 与 初始 压力 的 比值 P,/P; 来 描述 ,以 黑色 到 白色 
分 布 的 渐变 灰 度 图 呈现 。 如 图 5(a) 所 示 , 当 细胞 经 
过 两 微 柱 间 时 ,细胞 左 侧 受 到 较 大 压力 ,暂时 产生 了 
阻塞 。 同 时 ,由 于 细胞 膜 与 微 柱 间 液 膜 的 存在 ,接触 
处 细胞 膜 流速 没有 完全 为 零 , 微 柱 的 黏 清 效果 不 显 
著 。 图 5(b) 显 示 了 细胞 经 过 区 域 下 时 ,细胞 在 流体 
的 囊 挟 下 与 微 柱 碰撞 ,并且 开 始 横向 迁移 。 此 时 ,由 
于 惯性 和 细胞 膜 与 内 液 的 刚体 差 ,细胞 内 部 左 侧 出 
现 低压 区 ,而 右 侧 则 因 细 胞 内 液 和 微 柱 的 挤 压 作 用 ， 
出 现 高 压 , 液 膜 变 薄 , 微 柱 黏 滞 效 果 增 强 , 细胞膜 流 
速 接近 为 零 ,细胞 剪 切 率 增 大 。 在 图 5(c) 中 ,细胞 
继续 横向 迁移 ,此 时 细胞 的 流向 与 主流 道 方向 一 致 ， 
细胞 与 微 柱 接触 , 受 其 径 向 挤 压 力 和 黏 灌 力 的 影响 ， 
在 接触 处 产生 高 压 , 细 胞 膜 流速 也 变 组 ,细胞 剪 切 率 
较 大 。 


注 :细胞 直径 为 27.6 nm。 
图 5 细胞 运动 过 程 的 流 场 压力 分 布 
Fig.5 Flow field and pressure distribution during cell movement 


综 上 可 知 ,在 细胞 分 离 过 程 中 ,细胞 的 运行 轨迹 


取决 于 自身 大 十 和 DLD 芯片 几何 结构 ,而 结构 变化 
则 产生 于 其 与 微 柱 的 接触 过 程 中 。 因 微 柱 作用 力 以 
及 流体 压力 的 作用 ,细胞 膜 接触 部 位 的 结构 产生 显 
著 变 化 ,细胞 膜 与 微 柱 之 间 存 在 流体 液 膜 , 其 厚度 受 
压力 影响 ,从 而 影响 微 柱 对 细胞 的 黏 滞 和 剪 切 效果 。 


2.2 ”细胞 应 力 与 变形 分 析 


为 分 析 细 胞 膜 的 运动 特性 ,除了 要 了 解 细 胞 运 
动 过 程 中 的 结构 变化 \ 流 场 分 布 和 压力 情况 ,还 需 知 
道 细胞 膜 与 微 柱 接触 时 的 受 力 形变 情况 。 图 6 展示 
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了 细胞 运动 至 各 个 区 域 时 细胞 膜 所 受 应 力 和 相应 的 
变形 状况 。 其 应 力 为 依据 冯 “' 米 塞 斯 ( Von Mises) 
准则 ( 式 9) 求 得 的 等 效应 力 , 即 


T, = 了 Va -0,) +(0,-0) *(o,-o) ) 


(9) 

考虑 到 本 研究 所 使 用 分 子 动力 学 模型 为 周期 性 
边界 ,为 便于 展示 , 先 将 亚 区 域 的 两 片 细胞 膜 模型 按 
照 周期 性 合成 的 方法 ,整合 为 一 个 完整 的 细胞 模型 
(图 6) 。 同 时 ,依据 细胞 膜 的 模型 尺寸 (本 模型 细胞 
膜 宽度 为 3.4 nm) 和 细胞 运动 至 各 位 置 时 的 细胞 膜 
内 各 粒子 的 受 力 数 据 ,定量 划分 细胞 膜 的 受 力 宽度 
(黑色 部 分 为 细胞 膜 受 力 最 大 的 区 域 ) ,以 便 直观 展 
现 细胞 膜 的 受 力 情况 。 在 区域 ,只 有 细胞 膜 的 上 
下 边界 与 微 柱 发 生 接 触 , 因 而 细胞 膜 受 力 宽度 较 小 ， 
上 、 下 分 别 为 1.5、1.3 nm。 在 工区 域 ,细胞 右 侧 与 微 
柱 正面 碰撞 接触 ,细胞 膜 产生 较 大 变形 , 受 力也 相应 
增 大 , 受 力 宽度 为 2.6 nm, 而 细胞 膜 的 内 部 边界 ， 
碰撞 及 惯性 作用 和 刚度 差 , 细 胞 内 液 与 细胞 内 膜 分 
离 , 从 而 细胞 膜 内 边界 受 力 最 小 。 在 焉 区 域 ,细胞 向 
沿 着 主流 道 方向 运动 ,细胞 右 下 部 按 一 定 运 动 角度 
与 微 柱 发 生 接触 ,细胞 膜 变 形 和 受 力 都 很 大 , 受 力 宽 
度 为 3.3 nm, 超 过 了 1 了、 区域。 这 是 由 于 分 子 动力 
学 模拟 时 ,为 驱动 流体 流动 ,施加 了 外 势 场 ,在 DLD 
流 道 内 形成 了 压强 梯度 ,使 得 流体 流速 随 着 时 间 组 
慢 加 快 , 从 而 流体 压力 和 剪 切 应 力也 逐渐 提升 。 仿 
真 发 现 , 当 细 胞 直径 大 于 主流 道 宽度 时 ,细胞 膜 中 将 
产生 较 大 的 应 力 , 可 以 达到 GPa 级 别 。 以 往 的 文献 
表明 天” , 当 细 胞 膜 中 的 应 力 超过 0.1 GPa 时 ,细胞 
膜 有 可 能 产生 破裂 。 因 此 , 当 用 DLD 分 离 超过 主流 
道 尺寸 的 细胞 时 ,很 可 能 对 细胞 活性 产生 影响 。 

在 分 离 过 程 中 ,细胞 膜 在 微 流 道 的 整体 受 力 对 
比如 图 7 所 示 。 该 图 表示 了 归 一 化 后 的 细胞 膜 受 
情况 ,主要 对 比 了 细胞 膜 在 微 流 道 的 前 、 中 .后 3 个 
位 置 (对 应 图 4 的 工 .IT 、 亚 区 域 ) 和 不 同 微 柱 直径 时 
的 整体 受 力 情况 。 由 图 可 知 ,在 相同 微 柱 直径 时 , 细 
胞 膜 的 整体 受 力 在 微 流 道 的 前 .中 后 3 位 置 逐渐 增 
加 ;在 相同 的 位 置 时 , 微 柱 直径 为 14 nm 时 ,细胞 膜 
整体 受 力 最 大 ,为 20 nm 时 ,细胞 膜 的 整体 受 力 最 
小 。 结 合 图 6 的 模拟 结果 ,充分 考虑 流体 流速 增加 
对 细胞 膜 受 力 的 影响 ,再 一 次 表明 ,细胞 分 离 过 程 
中 ,细胞膜 受 力 的 主要 部 分 在 与 微 柱 的 接触 区 域 , 细 
胞 膜 变形 和 破坏 最 剧烈 的 部 分 也 在 该 区 域 。 
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注 : 细 胞 直径 为 27.6 nm。 
图 6 细胞 分 离 过 程 中 的 受 力 情 况 


Fig.6 Stress during cell separation 


细胞 平均 受 力 比 率 


14 16 18 20 22 
微 柱 直 径 /nm 


图 7 细胞 在 各 个 位 置 的 受 力 比率 


Fig.7 The stress ratio of cells at each location 


针对 DLD 装置 中 柔性 细胞 仿真 研究 的 难题 ,本 
研究 结合 粗 粒度 分 子 动力 学 模型 和 雷诺 相似 性 原 
理 , 尝 试 分 析 微 观 结构 对 细胞 变形 和 分 离 过 程 的 影 
响 。 选 取 粗 粒 化 分 子 模型 ,构建 了 细胞 DLD 分 离 的 
仿真 模型 ,分析 了 各 个 阶段 和 不 同 微 柱 直径 情况 下 ， 
细胞 在 DLD 芯片 中 的 运动 . 受 力 和 变形 机 制 。 

仿真 结果 发 现 , 在 纳 观 尺 度 下 的 DLD 细胞 分 离 
路 线 与 宏观 尺度 的 一 般 分 离 原理 相符 。 细 胞 的 分 离 
过 程 中 , 因 微 柱 作用 力 以 及 流体 压力 的 作用 ,细胞 膜 
接触 部 位 的 结构 产生 显著 变化 ,而 且 该 处 所 受 应 力 
最 大 。 细 胞 膜 与 微 柱 之 间 存 在 流体 液 膜 , 其 厚度 受 
压力 影响 。 分 析 了 不 同位 置 处 ,细胞 的 受 力 情况 ,以 
及 微 柱 直径 对 细胞 应 力 的 影响 。 

本 研究 的 工作 ,揭示 了 和 柔性 颗粒 通过 DLD 装置 
的 运动 细节 ,为 DLD 方法 的 研究 和 芯片 的 优化 提供 
了 参考 。 
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